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ABSTRACT 

On 2007 October 29 the outbursting comet 17P/Holmes passed within 0.79" of a background star. 
We recorded the event using optical, narrowband photometry and detect a 3% to 4% dip in stellar 
brightness bracketing the time of closest approach to the comet nucleus. The detected dimming 
implies an optical depth r 0.04 at 1.5" from the nucleus and an optical depth towards the nucleus 
center r n < 13.3. At the time of our observations, the coma was optically thick only within p < 0.01" 
from the nucleus. By combining the measured extinction and the scattered light from the coma we 
estimate a dust red geometric albedo pd = 0.006 ± 0.002 at a = 16° phase angle. Our measurements 
place the most stringent constraints on the extinction optical depth of any cometary coma. 
Keywords: comets: individual (17P/Holmes) — methods: data analysis — techniques: photometric 
— opacity — occultations 



1. INTRODUCTION 

An extraordinary explosion in late 2007 brought comet 
17P/Holmes (hereafter, 17P) back to the limelight of sci- 
entific interest, more than one hund red years after it was 
discovered by Edwin iHolmesI (|1892[ ). The discovery and 
early i nterest were tr iggered by an earlier, double out- 
burst (|BarnardlfT896l) . Since then and until the 2007 
event, 17P had received little attention due to its rather 
ordinary dynamical and photometric properties. 17P is a 
Jupiter-family comet, with an orbit between Mars (per- 
ihelion q = 2.1 AU) and Jupiter (aphelion Q = 5.2 AU) 
that has now experienced two explosive events separated 
by a comparatively quiet period of over 100 years. Here 
we report data taken shortly after the 2007 outburst to 
constrain the properties of the dust coma of this comet. 

Like most active comets, 17P appeared optically bright 
owing to scattering from dust particles expelled from the 
nucleus by gas drag. In principle, measurements of ex- 
tinction from cometary dust particles can be made by 
measuring the brightness of field sta rs when projected be - 
hind the coma of a passing comet (jCombes et al.llT983f ). 
Measurements of cometary extinction are potentially im- 
portant as, when combined with measurements of the 
scattered light, they can provide direct estimates of the 
ensemble dust albedo. In extreme cases, extinction ef- 
fects in cometary comae might influence the radiation 
budget at the nucleus and so influence the mass loss rate 
caused by sublimation. 

Several such m ea surem ent s have been re ported 
(lLarson fc A'Hearnl 119841: iLecacheux et all 119841: 
lElliot et alj (l995: Fer nandez et al.lll999t) but the results 
are mostly negative or difficult to interp ret, as a result of 
uncert ainties in the photometric data. iFcrnandcz et al.1 
(|1999( ) reported an extinction of about 20% in bright 
comet C/Hale-Bopp but, unfortunately, observed 
through highly non-photometric skies that bring the 
significance of the dimming into question. Even deeper 



Table 1 

Journal of Observations of 17P. 



UT Date 

Time of closest approach 

Heliocentric distance, R 

Geocentric distance, A 

Phase angle, ct 

Rate of motion 

Weather 

Telescope 

Instrument 

Pixel scale 

Seeing 

Filters (Exp. Time) 



2007 October 29 
UT 12:28:42 
2.455 AU 
1.627 AU 
15.6° 

26.4" hr" 1 
Photometric 
UH 2.2 m 
TEK 

0.219"/pixel 
0.9 ±0.1" 

BC (120 s), RC (25 s), R (1 s) 



events were reported in 95P/Chiron by lElliot et"aH 
(|1995|) and interpreted by them as resulting from 
extinction in narrow, near nucleus dust jets. Extinction 
from the diffuse coma of Chiron was not detected. 

The bright coma of exploding comet 17P provided 
an opportunity to search for measurable extinction of 
light from field stars. The comet impulsively ejected (2- 
90)xl0 10 kg of dust into the coma, with a peak pro- 
duction rate 3. 5 xlO 5 kg/s on UT 2007 Oct 24.54±0.01 
(|Li et al.l 120 111 ). These authors inferred, from sur- 
face photometry alone, an optically thick region near 
the nucleus subtending 0.1 arcseconds when the comet 
was at peak brightn ess. Extinction was reported by 
iMontalto et al.l (|2008f) four days later and at much larger 
angular distances from the nucleus (25 to 180 arcsec), 
where surface photometry indicates i mmeasurably lo w 
line of sight optical depths <2.5xl0" 3 (|Li et al.ll20TlT ). 

2. OBSERVATIONS 

Observations of 17P were obtained using the University 
of Hawaii 2.2-m telescope situated atop Mauna Kea. The 
telescope was equipped with a Tektronix TEK charged 
coupled device (CCD), which holds 2048x2048 pixels 
each 0.219" x 0.219" on the sky, and is read out in ~35 
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Figure 1. RC-band snapshots of the stellar appulse of 17P and star S. Differential photometry between stars S and C revealed dimming 
of the former as it passed behind the central coma of 17P. These are sub-frame sections of the TEK 2k CCD. 



s. Our measurements were obtained through the Halc- 
Bopp narrowband blue continuum (BC; A = 4453 A, 
AA = 61 A FWHM) and red cont inuum (RC; A = 7133 
A, AA = 5 8 A FWHM) filters (jFarnham et all l2000t 
Uewittl [20041) . For absolute calibration purposes we ob- 
tained a number of broadband R (Kron-Cousins-type fil- 
ter) images of 17P and Landolt standard stars through- 
out the night. To avoid saturation of the bright central 
region we were forced to use an integration time of 1 s, 
shorter than the minimum advised 5 s for the TEK cam- 
era. Integrations shorter than ~5 s suffer from a system- 
atic spatial pattern due to the finite shutter time that is 
more noticeable towards the edges of the CCD. Using 1 
s dome flatficld images we found that within our region 
of interest (near the center of the CCD as we tried to 
image the full circular coma) the shutter pattern uncer- 
tainty amounts to < 1%. The data were calibrated using 
bias frames and flatficld images obtained from dithered, 
median-combined images of the twilight sky. The night 
was photometric with very stable seeing, 0.9 ± 0.1". We 
set the telescope to track the non-sidereal rate (26" hr _1 ) 
of 17P and, as a result, stellar trailing in the longer BC 
exposures is comparable to the seeing. A journal of the 
observations can be found in Table [1] 



3. PHOTOMETRY 

We identified several field stars in the projected path of 
the comet and selected the one (TYC 3334-166-1) having 
the highest brightness (R = 11.0 ± 0.1 mag) and mini- 
mum impact parameter (p — 0.79") for this study; we 
refer to this as star S. Figure [1] shows some representa- 
tive images in RC band. By 2007 October 29, the coma 
of 17P had expanded to fill the entire 7 .5' field of view 
of the TEK detector (jHsieh et al.ll2010D . Consequently, 
each CCD pixel included contributions from the back- 
ground sky and from the 17P coma. To measure small 
changes in the flux from star S it was necessary to sub- 
tract the contribution from the sky and the coma of 17P. 
We attempted this subtraction in two different ways. 

Firstly, we employed pair-wise subtraction of near- 
consecutive frames aligned on the comet nucleus. Before 
subtraction, the frames were scaled by dividing each by 
the median pixel value measured in a 0.7' x 0.7' square 
region centered on the comet nucleus. The point spread 
function (PSF) of the star overlapped slightly in consec- 
utive frames (17P moved ~3.9" between frames) so we 
used every second frame for pair subtraction. A second 
method consisted of constructing "starless" frames by 
median-combining sets of three near-consecutive expo- 
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Table 2 

Model Optical Depths for 17P. 
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Figure 2. Lightcurve of star S versus projected radius from 
the nucleus of 17P through filters BC (blue circles) and RC (red 
squares). Negative projected radii indicate approach. The flux of 
star S was measured relative to that of star C. Overplotted as a 
solid black line is the inverted coma brightness profile, vertically 
shifted and scaled to fit the lightcurve of star S. 




Projected Radius [arcsec] 

Figure 3. Surface brightness of the coma of 17P versus projected 
distance from the nucleus along the path of star S. Overplotted 
as dashed lines are slopes proportional to p m , for m = —0.35 and 
m = — 1. The inset shows the slope of the profile as a function of 
projected distance. 
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Note. — Columns are (1) model for n(r) where "from SB" 
indicates that n(r) is inferred from coma surface brightness profile, 
(2) assumed q if constant, (3) assumed q for coma inner to 1", not 
resolved by our photometry, (4) measured optical depth at 1.5" 
impact parameter, (6) model optical depth towards nucleus. 



sures. Each median frame was then subtracted from the 
central of the medianed exposures. As before, the median 
of three consecutive exposures retained artifacts due to 
the star, so we used sequences 1-3-5 to construct the me- 
dian to be subtracted from frame 3. Medians of 5 frames 
(median of frames 1-3-5-7-9 subtracted from frame 5) 
produced similar results to 3-frame medians. Neither of 
the methods was ideal due to the fast-expanding coma, 
to slight changes in the seeing, and to imperfect align- 
ment of the nucleus in different frames. Both methods 
left traces of the coma that were strongest near the nu- 
cleus, but that amounted to at most, 1% of the flux from 
star S. The following results are based on the pair- wise 
subtraction method which we found to produce slightly 
better (less noisy) results. 

We performed circular aperture photometry centered 
on star S in each of the sky- and coma-subtracted images. 
Tests with increasing apertures indicated that an aper- 
ture radius of 3.3" produced the most stable photometry 
when star S projected far from the nucleus. Because of 
the subtraction, the flux surrounding the star should be 
zero, with a combined uncertainty that is approximately 
Gaussian. However, we chose to subtract the sky mea- 
sured in an annulus surrounding the star to remove local 
residual subtraction offset. To estimate the uncertainty 
in the flux from star S at position i we used the standard 
deviation of the fluxes measured at that same position i 
in the N — 1 frames (N = 15) for which the star is not 
present in the aperture. In addition, to confer protection 
against small seeing variations and other time-dependent 
variations we employed relative photometry. For this, we 
used another star (C, see Fig.HJ that lies further (>2.3') 
from the nucleus, which we measured in the same way 
as star S. Star C (R ~ 13 mag) lies distant enough from 
the bright central region of the coma that it should suffer 
much smaller extinction than star S (|Li et al.ll2011[ ). We 
computed the ratio of the fluxes of stars S and C as a 
function of time. 

Figure [2] shows the relative flux of stars S and C versus 
projected radius from the nucleus through the BC and 
RC filters. We divided the curve by its median value 
so that it explicitely represents the fractional dimming 
of star S. Near closest approach to the nucleus of 17P, 
the star dims by about 3% in BC (at 1.7" from the nu- 
cleus) and about 4% in RC (at 1.5"). The nearest BC 
and RC measurements occur on opposite sides of the nu- 
cleus and, for both bands, the brightness dip is visible in 
two separate measurements bracketing the nearest pro- 
jected distance from 17P. The two BC points are 0.6a 
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and 1.7cr below the median, while the two RC points are 
5.1cr and 4.0(7 below the median. Assuming that each 
pair of measurements is uncorrelated (i.e., their proba- 
bilities can be multiplied) the combined significance of 
the two dips is 2. Oct in BC and 6.6ct in RC. We note that 
both bands show a consistent trend of stronger dimming 
on the post-appulse side of the nucleus. This is true even 
though for RC the closest post-appulse measurement lies 
slightly farther from the nucleus (at ~ 2.5") than the 
nearest pre-appulse point. 

The observed dimming can be converted to an optical 
depth r using 



r(p) = - HF(p)/F ] 



(1) 



where p is the impact parameter at which the dimming 
is detected, F{p) is the star S flux at p and Fq is the 
unimpeded star S flux. Substituting 0.97 and 0.96 into 
equation (p} we obtain optical depths r(1.7") = 0.03 in 
BC and r(1.5") = 0.04 in RC. In the remainder of the 
paper we use the higher significance RC optical depth. 

4. DISCUSSION 

4.1. The optical depth to the nucleus 

Assuming a spherically symmetric coma composed of 
dust particles of radius a with a radial number density 
dependence n(r), the optical depth along a line-of-sight 
to the center of the nucleus, r„, can be related to the 
optical depth at a projected distance p by 



(l/2)J™na 2 n(r)dr 



t (p) Jo t™ 2 n(r')dl 



(2) 



where r' = \Jl 2 + p 2 , I is measured along the line-of-sight 
to projected radius p, r n is the radius of the nucleus, and 
the (1/2) factor accounts for the opaqueness of the nu- 
cleus. We can use the observed coma surface brightness 
profile (surface brightness versus projected radius from 
the nucleus, p) to infer the dust number density profile, 
n(r). If we assume the latter varies proportionally to r q 
then by integration the former will vary proportionally 
to p m \ where m = q + 1. Figure |3] shows the near-nucleus 
surface brightness profile from our RC data measured on 
a 3"-wide band centered on the path of star S across the 
coma. The profile is steepest near the nucleus, where it 
approaches the canonical slope p^ 1 of a spherically sym- 
metric, steady state coma, and is on average cx p~ - 35 
in the inner 30" of the coma. Using n(r) oc r q with the 
maximum slope q = — 2 (from m = —1), and repl acing 
r (1.5") = 0.04 and r n = 1.7 km (jLamv et al.ll2009t) into 
equation ([2]) we obtain r„ = 13.3. If, instead, we use the 
average inner coma slope q = —1.35 we obtain t„ = 0.19 
(Table [2]). The former, larger estimate can be regarded 
as an upper limit to the optical depth towards the nu- 
cleus. Steeper slopes than m = — 1 are generally due to 
solar radiation pressure which has insuff icient time to act 
in the inner 1" (~1000 km) of the coma (|Jewitt &: Meechl 
[1981 . 

Shown as an inset in Fig. [3] is the measured surface 
brightness profile slope as a function of projected radius 
p, which varies from m ~ —0.8 at p < 3" to m ~ —0.3 
at p ~ 1 2". These slopes ma t ch the m ~ —0.3 to —0.2 
found by iStevenson fc Jewittl (|2012ft at distances 10000 
to 20000 km (8" to 16" in our data) from the nucleus in 



data obtained a few days to a few weeks after our own ob- 
servations. We can use the measured radial dependence 
of m to solve numerically the integrals in equation ([2]) 
and improve our estimate of the optical depth to the nu- 
cleus, r n . Table [2] summarizes our results where, again, 
we used r(1.5") = 0.04 and r„ = 1.7 km. Because we 
only measure m to within approximately 1" of the nu- 
cleus we need to assume the value of q in the innermost 1" 
of the coma. We consider two possibilities: the limiting 
case q(p < 1") = —2, and a constant q(p < 1") = q{l"). 
In both cases we obtain optical depths to the nucleus 
only slightly larger than unity, falling off to t = 1 at 
p ~ 0.01". 

Our measur ement is consistent with the findings of 
iLi et all (|201lD who used coarser spatial resolution pho- 
tometry on data taken around 2012 October 24.5 to 
conclude that the coma of 17P was optically thick only 
within a tiny central region, ~ 0.1". They find a peak 
optical depth towards the nucleus r n ~ 65 on October 
24.8. Another dete ction of extinct i on by the coma of 
17P was reported bv lMontalto et all (|2008l ). using an en- 
semble of 20 stars located between 25" and 180" from 
the nucleus of the comet. They infer optical depths 
0.3 < r < 0.5, two to three orders of magnitude larger 
than the 1 x 10 -4 < r < 4 x 10 -3 implied by our mea- 
surements in the region 25" < p < 180". 

4.2. The coma dust albedo 

The extinction obtained in §4.11 can be combined with 
the measured scattered light, to provide a direct estimate 
of the coma dust albedo. We here use the R band images 
mentioned in Assuming that the inner coma of 17P is 
spherically symmetric then the optical depth at projected 
radius p from the nucleus can be written 



r(p) 



A, 



7T (2 p 8p) A 2 



(3) 



where ir(2p5p) A 2 is the projected area of the annulus 
centered on the nucleus with inner radius p — Sp/2 and 
outer radius p + 5p/2, and Ad is the total effective scat- 
tering cross-section of the dust grai ns within the a nnulus. 
The latter can be calculated from (jRussellll 19161) 



Pd 



A d = 2.24 x io^ttIO - 41 "^-™^ 1 



(4) 



where pd is the red albedo of the dust grains, tuq = 
—27.11 is the apparent red magnitude of the Sun, and 
rrid(p) is the absolute red magnitude^ of the dust con- 
tained within the annulus defined above. Equations ([3]) 
and (2|) can be combined to yield 



Pd 



2.24 x 10 22 1 0°- 4 [ m o~ md ( p )J 
2p5pA 2 T(p) ' 



(5) 



We take p = 1.5", Sp/2 = 0.219" (1 pixel), t(1.5") = 0.04 
and calculate the magnitude of the coma in the annulus 
using 



m d (1.5") = -2.51og 10 (l0 
= 7.17 ±0.47 mag 



0.4m out 



10" 



(G) 



The absolute magnitude is the apparent magnitude reduced 
to unit heliocentric and geocentric distances and to zero degrees 
phase angle. 
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where m out = 6.04 ± 0.11 mag and TO; n = 6.51 ± 0.11 
mag are the absolute red magnitudes measured within 
the outer and inner radii of the annulus. Replacing those 
quantities into equation [5] we obtain pd = (6 ±2) x 10~ 3 . 

Our estimated dust albedo is significantly lower than 
an earlier, i ndependent measurem ent using a different 
technique by llshiguro et al.l (|2010l ). Those authors used 
combined optical and infrared observations of 17P taken 
on 2010 October 25-28 to find an albedo 0.03 < p d < 
0.12. Both estimates were obtained at sim i lar ph ase an- 
gle, a = 16°. Interestingly, llshiguro et afl (|2010l ) report 
a decrease in albedo with time, from pd = 0.12 ± 0.04 
on October 25 to p d = 0.03 ± 0.01 on October 28. Our 
estimate follows the trend but we find a 5 x lower albedo 
just a day later, on October 29. 

Earlier attempts to estimate the reflectivity of 
comctary dust grains using stellar appulses have also 
found surprisingly low albedos. lLarson fc A'Hearnl 
(|1984|) observed comet Bowell (1980b) and found a dust 
albedo of pd ~ 1.5 x 10 -3 . They suggested that pho- 
tons might go through two reflections (0.0015 « 0.04 2 ) 
before escaping the dust grain to reconcile their find- 
ings with the typical 4% albedo of cometary surfaces. 
iFernandez et al.l (|1999|) reported an albedo pd ~ 0.01 for 
the coma dust of comet Hale Bopp. 

Such low albedos are not unprecedented. Spacecraft 
imaging of comet Borrelly revealed very dark regions 
with geomet ric albedos p < 0.01 at wavelengths 0.5 p,m < 
A < 1.0 /im (|Soderblom et al.l 12001: iNelson et al.ll200"l . 
High-porosity carbon-based aerogels are also known to 
have very low albe dos over a broad range of wavelengths 
(Merzbachg3 120011 ) . The least reflective carbon aero- 
gels (reflectivities as low as 0.003) arc made of ~1 ran 
amorphous carbon particles and ~10 ran voids resulting 
in bulk porosities larger than 90%. Indeed, theoretical 
studies predict that the albedo of porous grain aggre- 
gates of small particles should decrease with increasing 
bulk porosity (IHage fc Greenberel 119901) . reaching sub- 
1% reflectivities at void volume fraction of about 70% 
(|Napier et al.ll200"4T ). The low reflectivity of carbon aero- 
gels is due in part to strong absorption from amorphous 
carbon but results mostly from the rough, highly porous 
physical structure of the material (silica aerogels are also 
highly porous but are smoother on 1 pm scales). Our 
observations and those of others may be suggesting that 
comet dust is highly porous and, in the case of 17P, pos- 
sibly carbon-rich. We note that a significant amount of 
amorphous carbon is required to explain Spitzer spectra 
of 17P taken on 2 007 November 10, only a few days after 
our observations (jReach et al.H2010l) . 

5. CONCLUSIONS 

We used narrowband photometry of a background star 
that passed within < 1" of the nucleus of out bursting 
comet 17P/Holmes to constrain the extinction optical 
depth of the cometary dust coma. Our conclusions are 
as follows: 

• We detect a 3 to 4% dimming of the background 



star depending on wavelength that coincides with 
the nearest projected distance to the cometary nu- 
cleus. If caused by extinction due to coma dust, 
this dimming implies an optical depth r = 0.03 to 
0.04 at projected radius 1.5" from the comet nu- 
cleus. We infer that the coma was optically thick 
only in a tiny region, ~ 0.01" in radius, surround- 
ing the nucleu s. Our m e asure ments are compatible 
with those bv lLi et all (120111) but d i fficult to rec- 
oncile with those bv iMontalto et all (|2008| ). 

• We combine the measured extinction with pho- 
tometry of the coma to estimate a dust albedo of 
(0.6 ± 0.2)%. Our albedo estimate is 5x lower 
than an independent measurement using optical 
and t hermal observations of light scattered by the 
dust (jlshiguro et al.ll2""010D , but it is consistent with 
the darkest regions on comet Borrelly, and with the 
dust being composed of very small, highly porous 
grains of carbonaceous composition. 
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